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ABSTRACT 


Is  there  an  upper  limit  to  the  number  of  elements  or,  equivalently,  samples,  that  can  be  usefully 
employed  for  direction  finding  or  nulling  a  set  of  incoming  wavefronts  within  a  fixed  size  aperture? 
Mathematical  theory  states  that  this  number  isunlimited,  but  this  result  requires  a  perfect  knowledge  of 
the  wavefront  correlation  matrix  (i.e.,  a  noise-  and  error-free  environment)  and  infinite  precision.  To 
answer  the  question  realistically,  one  must  consider  the  element-to-element  correlation  matrix  resulting 
from  a  set  of  incoming  wavefronts  impinging  on  an  array  of  antenna  elements.  This  matrix  becomes  ill- 
conditioned  as  the  number  of  elements  within  the  aperture  increases. 

This  report  shows  that  each  eigenvalue  associated  with  the  correlation  matrix  equals  the  sum  of  the 
powers  of  all  the  wavefronts  projected  onto  the  dimension  represented  by  its  corresponding  eigenvector. 
If  the  eigenvalue  is  of  negligible  power  and  easily  lost  in  noise  and  other  causes  of  data  corruption,  then 
the  essential  information  associated  with  that  eigenvector,  and  the  dimension  it  represents,  is  lost. 
Elements  corresponding  to  the  number  of  small  eigenvalues  are  of  little  value  for  the  given  set  of 
wavefronts.  The  results  must  be  generalized  by  choosing  a  set  of  wavefronts  that  encompasses  the  entire 
range  of  interest  so  that  if  an  element  is  of  questionable  value  for  this  “testing  wavefront,”  then  it  has 
doubtfiil  value  for  any  set  of  wavefronts.  Using  the  testing  wavefront,  an  upper  bound  can  be  determined 
to  the  number  of  useful  elements  that  can  be  employed  within  a  specified  fixed  aperture. 
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1.  INTRODUCTION 


How  many  antenna  array  elements  (samples)  are  beneficial  for  direction  finding  or  nulling  with  a 
fixed  aperture  size?  Mathematically,  additional  elements  are  beneficial;  however,  perfect  knowledge  of 
the  ideal  correlation  matrix  and  infinite  precision  are  required.  With  finite  precision,  finite  time  averaging, 
finite  signal  to  noise,  finite  dynamic  range,  radio  firequency  interference,  antenna  element  coupling, 
nonuniform  directional  background  noise,  and  pattern  distortions,  the  number  of  useful  elements  is 
compromised.  To  answer  the  question  realistically,  the  correlation  matrix  resulting  from  a  set  of  incoming 
wavefronts  is  considered;  the  matrix  becomes  ill-conditioned  as  the  number  of  elements  increases  within 
the  aperture.  Each  eigenvalue  associated  with  the  correlation  matrix  corresponds  to  the  power  sum  of  all 
the  wavefronts  projected  onto  the  dimension  represented  by  its  corresponding  eigenvector.  If  the 
eigenvalue  is  of  negligible  power  and  easily  lost  in  noise  and  other  causes  of  data  corruption,  then  the 
essential  information  associated  with  that  dimension  is  lost.  By  reducing  the  number  of  elements  to 
correspond  to  the  number  of  significant  eigenvalues  and  suitably  positioning  the  remaining  elements 
within  the  specified  aperture,  little  information  regarding  the  set  of  wavefronts  is  lost. 

The  result  is  generalized  by  postulating  a  set  of  wavefronts  that  encompasses  the  entire  range  of 
interest;  if  an  element  has  questionable  values  for  this  wavefront  set,  designated  the  testing  wavefront, 
then  it  has  doubtful  value  for  any  set  of  wavefronts.  Using  the  testing  wavefront,  an  upper  bound  can  be 
determined  to  the  number  of  useful  elements  that  can  be  employed  within  a  specified  fixed  aperture. 

Linear  arrays  of  fixed  aperture  are  considered,  then  generalized  to  two-dimensional  circular  arrays. 
This  report  is  confined  to  linear  and  circular  arrays,  but  the  technique  can  be  applied  to  more  general  (e.g., 
three-dimensional)  arrays. 
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2.  BACKGROUND 


Considered  here  is  a  linear  array  of  M  isotropic  antenna  elements  positioned  at  locations  dm  along 
the  x-axis.  The  M  vector  of  baseband  receiver  outputs  at  a  sampling  time,  t,,  due  to  the  kth  narrowband 
received  signal  of  wavelength  A  and  planewave  direction  ^  off  broadside  to  the  axis,  is 


where  the  complex  scalar  represents  the  amplitude  and  phase  of  the  received  signal  at  element  one 
(using  reference  element  w  =  1  so  that  di  =  0).  is  the  M-dimensional  vector  of  samples  that  is 
referred  to  as  the  “response  vector”  and  defined  by 

-  (1)  ^kidi),  J k{d-i),.-;zl k(di^Y  , 


where,  in  turn, 


4(^4)  =  Qxp(-j2^d„  sin  dk/X) 


so  that  ^k(d\)  =  1.  The  Sk(t,)  are  modeled  as  zero-mean,  independent  (in  both  k  and  /)  random  sequences  of 
stationary  statistics. 

Uniform  linear  arrays  are  considered  first,  then  generalized  to  nonuniform  and  two-dimensional 
arrays.  For  the  uniform  array,  dm  —  (rn—l)d,  where  d  is  the  common  separation  between  elements.  In  this 
case,  z! kidj)  becomes  z*  and  z^ *(<4)  =  z*"^',  where 

Zk  =  exp(-j2zrd  IX  sin  6^)  . 


The  response  vector,  l^(zi),  then  simplifies  to  V(zk),  where 


This  latter  vector  F(z*)  is  called  a  Vandermonde  vector.  In  this  case  £*  represents  F(zi).  Any  distinct  set  of 
M  or  fewer  Vandermonde  M-dimensional  vectors  is  linearly  independent.  In  order  that  no  two  off- 
broadside  angles.  Ok,  share  the  same  value  of  Zk,  d  must  be  no  larger  than  X2. 

If  there  were  K  signal  wavefronts  {K<  M)  impinging  on  the  array,  the  vector  of  receiver  outputs  at 
time  t\  would  be 
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(1) 


K{t,)=n{t,)+'t 

k=l 

where  «(//)  is  the  M  vector  of  internal  receiver  noise  and  other  noise  sources,  which  for  this  analysis  are 
assumed  independent  from  receiver  to  receiver  and  of  average  power  level  &  in  each  receiver. 

The  correlation  (covariance)  matrix,  R,  is  obtained  by  averaging  the  Af  x  M  outer  product 
where  the  H  operator  is  the  conjugate  transpose.  It  is  assumed  that  the  signal  wavefronts  are 
all  uncorrelated  so  that 

=  0  iork^j  .  (2) 

With  this  signal  model,  the  expectation  for  the  R  matrix  becomes 

R  =  E{X(tmtf}  =  VPV^  +  dl  , 

where  P  is  a  AT  x  AT  diagonal  matrix  with  diagonal  terms  Pk,  which,  in  turn,  are  given  by 

=  forA:  =  i,2,...,is:  .  (3) 

F  is  an  M  X  AT  matrix  having  columns  that  are  K  Vandermonde  vectors,  each  column  corresponding  to  a 
signal  direction.  The  signal-dependent  component  of  the  R  matrix,  S,  can  be  spectrally  decomposed  to 

S  =  VPV''  =  'Z  EkKkEf  (4) 

k=l 

to  obtain 

5  =  FPF"  =  X^*  •  (5) 

k=I 

Note  that  the  subscripted  As  refer  to  eigenvalues,  whereas  the  unsubscripted  Prefers  to  wavelength.  The  K 
Vandermonde  vectors  span  a  space,  Sk,  of  dimension  K.  The  M  eigenvectors,  span  a  space  Sm  of 
dimension  M.  The  space  spanned  by  Sk  for  AT  <  Af  is  a  subspace  of  Sm-  The  space  Sk  is  referred  to  as  the 
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“signal  subspace”  that  is  identical  to  the  space  spanned  by  the  K  eigenvectors  of  S  having  nonzero 
eigenvalues.  The  remaining  M-K  dimensional  subspace  is  the  noise  space.  Since  all  the  eigenvectors  of  a 
Hermitian  matrix  are  mutually  orthogonal,  the  Vandermonde  vectors  are  each  orthogonal  to  the  noise 
subspace.  The  noise  space  can  be  used  to  direction  find  up  to  M-1  wavefronts. 

In  an  ideal  situation,  S  is  obtained  fi’om  R  by  subtracting  C5^7.  In  reality,  there  is  no  R,  but  it  can  be 
estimated  with  a  sampled  version,  r!,  by  averaging  over,  say  L  samples, 

,  (6) 

L 


where,  for  the  error  matrix,  is 


E  = 


2 

Oi  £12  £im 
2 

€21  (y 2  •••  siM 

2 

Smj  £m2  •••  (7m 


(7) 


and  where  L  is  as  large  as  practical. 

For  the  ideal  case,  E  is  simplified  to  dl.  If  it  is  assumed  that  the  K  wavefronts  are  incoherent,  then 
the  signal  cross-terms  in  ^  average  down  to  near  zero,  leaving,  essentially,  the  signal  self-terms.  The 
independence  of  the  noise  terms  causes  both  the  noise  cross-terms  and  the  signal-noise  cross-terms  to 
average  down  to  nearly  zero  so  that  the  £mn  terms  should  all  be  close  to  zero. 

If  all  the  element  patterns  and  receivers  are  identical,  then  the  Cm  are  about  equal,  giving  R~r!.To 
ascertain  the  estimate  of  S,  a  value  must  be  determined  to  be  used  for  and  I  must  be  subtracted  fi'om 
1^.  In  the  ideal  case  this  step  is  relatively  easy  because  with  K<  M,  the  noise  eigenvalues  of  R  are  all  (} . 
Instead,  there  is  /?^,  and  since  E  has  nonzero  off-diagonal  terms  and  the  Om  are  not  identical  on  the 
diagonal,  the  choices  for  K  and  &  are  usually  not  obvious  [1-5].  With  one  nonzero  pair,  Bmn-  ^ 
each  eigenvalue  is  shifted  by  an  amount  proportional  to  e,„„  with  the  proportionality  factors  summing  to 
one  [6].  Critical  information  in  the  eigenvalues  and  eigenvectors  is  thereby  corrupted.  A  typical 
noise/distortion  floor  is  expected  to  be  in  the  range  of  30  dB  below  the  peak  signal-to-noise  ratio  (SNR). 
With  precision  calibration,  a  further  reduction  can  be  expected  on  the  order  of  15  to  20  dB,  which  is  the 
current  limit,  but  even  if  a  10  dB  improvement  is  allowed,  it  would  make  a  difference  of  no  more  than  one 
useful  sample.  As  a  result,  most  systems  will  have  difficulty  distinguishing  between  signal  and  noise 
eigenvalues  due  to  distortions  and  nonimiform  noise  when  the  are  about  40  dB  below  the  largest  signal 
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eigenvalue.  For  the  most  optimistic  situation,  another  10  to  15  dB  may  be  allowed,  which  might  admit 
one  or  two  more  useful  samples;  the  number  of  useful  elements,  or  samples,  is,  therefore,  limited. 
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3.  DETERMINING  SAMPLING  LIMITS 


Information  about  the  signal  waveforms  impinging  on  the  array  is  obtained  from  the  larger 
eigenvalues  and  the  corresponding  eigenvectors  of  the  S  matrix.  A  correlation  matrix  must  be  constructed 
that  ultimately  will  allow  the  determination  of  the  number  of  useful  samples;  i,e.,  the  useful  number  of 
elements  in  the  array  for  a  given  aperture  size  in  wavelengths.  A  Vandermonde  vector,  ii,  is  expanded  in 
terms  of  the  eigenvectors  of  S,  to  obtain 

M 

=  where  (8) 


SO  that  the  akm  are  the  projections  of  the  Ath  Vandermonde  on  the  mth  eigenvector.  By  definition, 

K=yLlsu^ 


so  that 


^  He  _  H 

— m  — m 


S/’.  5:*  t", 


k^l 


=  t.  Ptia"  E.)(rf  •£„)=i;Ak 


k=l 


k=l 


(9) 


Each  eigenvalue  is  equal  to  the  sum  of  the  projected  powers  of  all  the  Vandermondes  onto  its 
eigenvector  [see  Equation  (9)].  Small  eigenvalues  correspond  to  dimensions  in  which  there  is  almost  no 
power  from  any  wavefront.  In  an  ideal  case,  any  nonzero  level  is  satisfactory  for  extracting  information, 
which  is  not  true  in  the  real  environment,  where  S  is  not  known.  Note  that  if  eigenvalues  are  too  small, 
they  cannot  be  determined  with  good  accuracy,  and  the  coordinates  corresponding  to  their  eigenvectors 
provide  no  useful  information  [see  Equation  (9)].  Therefore,  fewer  coordinate  axes,  i.e.,  fewer  samples, 
would  be  as  useful,  less  complicated,  and  less  expensive. 

If  interest  is  in  a  specific  set  of  wavefronts,  the  corresponding  correlation  matrix  can  be  formed  and 
its  eigenvalue  rolloff  analyzed.  Generally,  however,  a  set  of  wavefronts  must  be  selected  that  generate  a 
matrix  S  that  allows  one  to  determine,  for  a  given  aperture,  coverage  region,  and  signal  spatial 
distribution,  the  maximum  number  of  useful  samples  that  can  be  taken  imder  some  total  noise/distortion 
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level  assumption.  The  resulting  eigenvalue  rolloff  then  gives  an  upper  bound  to  the  number  of  useful 
elements/samples  for  a  particular  signal  spatial  distribution. 

Given  that  it  is  not  desirable  to  artificially  limit  the  number  of  significant  eigenvalues  due  to  lack  of 
signals,  a  set  of  wavefronts  must  be  selected  with  K>  M  and,  preferably,  K  »  M,  as  well  as  choosing 
signal  wavefronts  that  will  not  bias  the  results.  To  achieve  this  end,  the  signal  directions  and  the  signal 
power  levels  are  made  uniform  (P  =  7)  so  that  identical  signal  plane  waves  uniformly  cover  the  entire 
region.  The  resulting  set  of  wavefronts  is  referred  to  as  the  “testing  wavefront.”  The  intent  is  to 
manufacture  a  covariance  having  as  many  high  eigenvalues  as  possible.  With  Mas  the  number  of  antenna 
elements,  A/ A  as  the  aperture  in  wavelengths,  and -W  to  W  (W<7if2)  as  the  coverage  region,  S  is 


■S=s  P,  K,  K“  , 


(10) 


k^/ 


with 


1 

2W 


(11) 


where  signals  are  separated  by  an  angle  S<p.  If  approaches  zero,  then,  since  KS(p-  W,  approaches 


1  ^ 


{m-n)A  . 
’2^-y - rrsin 


(12) 


The  Bessel  function  identity  [7,  #9.1.41  with  /  =  exp(^’ci)] 


^  V/Z) 


(13) 


can  be  used  in  Equation  (12)  to  obtain,  after  changing  the  order  of  integration  and  summation. 
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1 


oo 


s 


mn 


2W 


^  Itu  A{m-n)^ 

,  KM-\)  , 


W 

I 

-w 


dip 


2n  A{m-n)\  sin(p^) 

X{M-\)  J  pW 


(14) 


With  the  identity  [7,  #9.1.5] 


(15) 


Equation  (14)  can  be  written 


■Jo 


Ik 


A  (m  —  n) 


A 


+  2  E 

p^2 

p  even 


2k  — 

A 


(w  -  «)  I  sin  (p  W) 
(M-l) )  pW 


(16) 


and  note  that  for  W=  7i/2  (90°),  s„„  simplifies  to  the  Jo  term  only.  For  a  nonuniform  array,  with  the  wth 
element  located  at  d„  and  A  ^d^-du  the  result  in  Equation  (16)  is  generalized  to 


=  Jo 


2k 


+2  E  J, 


'lUfj  ^sinCf  r) 

/»  \U  m  Ct  n) 


p  =  2 

p  even 


V 


pW 


(17) 


Since  the  testing  wavefront  contains  fronts  from  more  directions  than  most  practical  scenarios, 
eigenvalues  are  expected  to  have  more  sources  contributing  to  their  value,  they  roll  off  more  slowly,  and 
an  upper  bound  to  the  count  is  being  produced. 

For  a  uniform  array,  s„„  =■  s^JM,  Al'k,  is  a  ftinction  of  three  parameters;  M  is  the  number  of 
samples  (or,  equivalently,  elements);  AIX  is  the  aperture  (in  wavelengths),  and  the  observation  range  is 
—W  to  W.  The  investigation  of  signal  space  behavior  and  its  eigenvalues  begins  for  the  cases  in  which  the 
array  is  uniform  and  W  -  7ti2  (90°)  and  generalizes  the  results  to  other  cases.  The  intent  is  that  the 
resulting  S  matrices  should  have  the  slowest  eigenvalue  drop-off  for  the  given  parameters,  hence  the  most 
measurable  eigenvalues.  Any  other  signal  arrangement  produces  a  set  of  eigenvalues  that  roll  off  more 
quickly. 
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4.  UNIFORM  ARRAYS 


Generated  first  are  the  S  matrices  for  W  =  /r  12(90°),  a  set  of  A/ As,  and  a  sequence  of  different  values 
of  M.  Eigenvalues  are  arranged  in  decreasing  order  for  each  case,  and  the  rolloffs  of  these  eigenvalues  are 
examined.  Figure  1  shows  A/1  =  1  and  M  varying  from  6  to  10.  Note  that,  as  far  as  they  can,  the  relative 
eigenvalues  follow  almost  identical  rolloff  sequences,  nearly  independent  of  M  Figures  2  and  3  show 
similar  results  for  A/1  =  2  and  A/ 1=3. 


344558-1 


Figure  1.  Eigenvalue  rolloff  versus  m  for  M  =  6-10,  A/1  =1,  uniform  array,  testing  wavefront. 


An  important  property  in  all  cases  is  that  as  M  increases  and  the  spacing  between  elements  falls 
well  below  AJ2,  the  S  matrix  becomes  increasingly  ill  conditioned;  with  increasing  M,  some  eigenvalues 
become  too  small  to  allow  an  accurate  calculation  of  the  required  projection  coefficients.  The  dimension 
number  beyond  which  the  rolloff  becomes  rapid  is  dependent  on  the  aperture  of  the  array  and  W;  the 
larger  the  aperture,  the  larger  the  dimension  number  at  the  knee  in  the  curve.  Since  the  knee  is  ever 
present,  a  definite  limit  exists  to  the  number  of  useful  elements  that  can  be  packed  into  a  fixed-size 
aperture. 

Holding  A/1  =  2  and  M=  11,  IF  is  reduced  to  5;r/12  (75°),  n:/3  (60°),  and  ;r/4  (45°)  and  results  are 
determined  with  the  change  in  W.  Figure  4  shows  that  the  drop-off  occurs  somewhat  earlier  as  W 
decreases,  but  the  results  are  not  particularly  sensitive  to  changes  in  W.  The  slowest  rolloff  occurs  when 
W  =  n/2  (90°),  i.e.,  the  maximum  number  of  useful  samples  for  a  given  threshold  level. 
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344558-2 


Figure  2.  Eigenvalue  rolloff  versus  m  for  M  =  7-11,  A/X  =2,  uniform  array,  testing  wavefront 


344558-3 


Figure  3.  Eigenvalue  rolloff  versus  m  for  M  =  9-13,  A/X  =3,  uniform  array,  testing  wavefront. 
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344558-4 


Figure  4.  Eigenvalue  rolloff  versus  m  for  M  ~  If  A/X  —2,  W  —  45^,  60°,  75°,  and  90°,  uniform  array,  testing 
wavefront 
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5.  NONUNIFORM  ARRAYS 


Each  element  location,  d^,  must  be  specified  for  the  nonuniform  array.  S  matrices  are  generated 
using  Equation  (17)  in  place  of  (15).  While  the  S  matrices  were  symmetric  and  Toeplitz  for  uniform 
arrays,  these  general  S  matrices  are  not  Toeplitz. 

It  was  shown  earlier  that  the  W  =  nil  (90°)  cases  produce  the  slowest  rolloff  and  that  the  results  are 
not  sensitive  to  the  choice  of  W.  Initially,  different  cases  with  W  —  it  IT.  (90°)  are  considered  for  the 
nonuniform  array.  To  guarantee  that  these  arrays  cover  the  full  aperture,  element  one  is  at  position  zero 
and  element  M  at  position  A.  Even  though  most  nonuniform  arrays  are  designed  to  thinly  or  sparsely 
sample  an  aperture,  it  is  being  oversampled  here  to  show  that  the  sampling  results  obtained  for  uniformly 
spaced  arrays  are  not  dependent  on  uniform  sampling.  Two  well-known  nonuniform  spacings  are 
considered — log-periodic  and  random. 

5.1  LOG-PERIODIC  SPACING 

Log-periodic  (more  precisely,  periodic  log-frequency)  spacing  specifies  most  element  locations  as 
scaled  by  a  constant  factor,  “a”  (a  >  1),  from  the  previous  location.  The  first  element,  du  is  at  zero.  The 
next,  d2,  is  arbitrary.  The  third  and  all  subsequent  elements  are  scaled  by  “d”  so  that  d-i  —  ad2,  d^  —  ad-i  — 
a^d2,  etc,  with  dM  =  a^'^d2  being  at  the  aperture  point  A .  The  separation,  d2-d2=(a-  \)d2  is  less  than  d2  if 
a  <2  and  would  then  be  the  smallest  spacing,  preferably  <  XI2. 

The  results  for  this  array  are  given  in  Figures  5-7.  The  same  behavior  as  the  uniform  array  is  seen 
except  that  the  rapid  fall-off  is  somewhat  greater  than  with  the  uniform  array,  indicating  that,  for  a  dense 
array  with  average  spacing  less  than  M2,  the  uniform  array  is  the  better  arrangement.  Figure  8  is  an 
example  of  the  eigenvalue  fall-off  for  different  Ws.  Results  are  very  similar  to  the  uniform  case.  Note  that 
log-period  spacing  is  usually  adopted  for  sparse  arrays  where  economy  in  the  number  of  elements  is  a 
primary  issue. 

5.2  RANDOM  SPACING 

Except  for  the  extreme,  elements  are  randomly  placed  within  the  aperture;  therefore,  if  there  are 
nine  elements,  two  are  fixed-end  positions,  with  seven  randomly  placed  between  them  for  a  uniform 
distribution.  Results  are  shown  in  Figures  9-11  and  are  almost  identical  to  those  of  the  log-periodic 
spacing  as  is  the  varying  W  example  given  in  Figure  12.  Both  nonuniform  arrays  exhibit  behavior  similar 
to  the  uniform  array.  There  is  clearly  a  practical  limit  to  the  number  of  useful  samples  for  linear  arrays. 
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0 


Figure  7.  Eigenvalue  rolloff  versus  mforM-  9-13,  A/X  =3,  log-periodic  array,  testing  wavefront 


344558-8 


Figure  8.  Eigenvalue  rolloff  versus  m  for  M  =  II,  A/X  -2,  W  ~  45^,  60^,  75^,  and  90°,  log-periodic  array,  testing 
wavefront. 
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344558-9 


Figure  9.  Eigenvalue  rolloff  versus  m  for  M  -  6-10,  A/X  =  1,  random  array,  testing  wavefront 


Figure  10.  Eigenvalue  rolloff  versus  m  for  M  =  7—11,  A/X  —  2,  random  array,  testing  wavefront 
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344558-11 


Figure  II.  Eigenvalue  rolloff  versus  mforM-  9-13,  A/X  -  3,  random  array,  testing  wavefront. 


344558-12 


Figure  12.  Eigenvalue  rolloff  versus  mfor  M  ^  II,  A/X  -2,  W  =  45°,  60°,  75°,  and  90°,  random  array,  testfr 
wavefront. 


6.  CIRCULAR  ARRAYS 


Arrays  with  two-  or  three-dimensional  apertures  can  also  be  investigated.  For  a  general  three- 
dimensional  array  with  elements  at  position  where 


and  n  is  the  direction  vector  (^is  the  spherical  azimuth  angle  and  ^the  spherical  elevation  angle) 


={  cos  ^  sin  sin  (p  sin  0,  cos  6)  , 


(18) 


the  response  vectors  become 


where,  in  turn, 


i^(4)=(i,4(^),4(^),...y«)^  , 


4(^L)  =  exp(-y2;r(ci„-n) ) 


Then  the  terms  of  S,  s^m  become 


1  1  ^ 

— L  f  [ 

iwuLl 


(19) 


where  0  to  f/is  the  range  on  the  elevation  angle  6-  ujjn,n),  Uy(m,n),  and  uJim,n)  are  defined  as  satisfying 
2  7C 

~r Him =  Uxi^> «) cos ^ sin ^  +  Uyijn, n) sin ^ sin ^  n) cos 6  .  (20) 

A 

In  particular,  a  radius  r  circular  array  in  the  x-y  plane,  with  M  uniformly  spaced  elements,  gives 

^^  =  r{cos{m-\)27t/M,sm.(m-\)l7t/M,0)  ,  (21) 

and  the  ujjn,ri)  and  Uy(m,n)  become 
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(22) 


Ux{^>  n)-2K~  [cos  27t{m  -  X)/M  -  cos  2K{n  -  1)/M] 
A 


and 


i 

Uyitn,  n)-27r  —  [sin  2K{m  -  \)/M  -  sin  27t{n  -  \)/M] 


(23) 


with  Ux{m,n)  =  0.  For  the  uniform  circular  array,  the  become 


1  1 

2WU 


u  w 


0  -w 


(24) 


For  a  numeric  generation  of  the  s„n  used  for  general  W  and  U, 
1  1  u/Sd  w/Sip 


^VV  U  k  =  -W/5<l> 

where,  specifically, 


(25) 


Ux(jn>n)^'^^{xm-x„)  and  Uy(rn,n)^^^{y^-yJ 


(26) 


Using  Equation  (25)  for  a  numerical  evaluation,  the  S  matrix  can  be  filled  and  the  behavior  of  its 
eigenvalues  examined. 

The  focus  now  turns  to  two  circular  arrays.  The  first  has  all  its  elements  uniformly  distributed  on 
the  circumference.  The  second  (referred  to  as  “circular-0,”)  is  similar  to  the  first  except  that  one  element 
is  placed  in  the  center  of  the  circle.  A  testing  wavefront  is  chosen  to  cover  a  ±Ji/2  (±90°)  wedge  in 
azimuth  and  either  a  7c/4  (45°)  or  %I2  (90°)  wedge  in  elevation.  Again,  the  testing  wavefront  is  expected  to 
produce  a  slow  eigenvalue  rolloff  and  an  upper  bound  to  the  number  of  useful  elements.  For  example,  the 
radius  of  the  circular  array  is  0.3X.  A  circle  with  a  0.3A,  radius  has  a  1.885X  circumference.  With  a  four- 
element  array,  the  elements  are  separated  on  the  circumference  by  approximately  'kl2.  As  such,  the 
eigenvalue  rolloff  is  expected  to  start  somewhere  between  five  and  six  elements.  The  resulting  rolloffs  for 
four-,  five-,  six-,  and  seven-element  arrays  are  plotted  in  Figure  13. 
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344558-13 


Circular  Array,  r  =  0,3% 


Figure  13,  Eigenvalue  rolloff  for  the  circular  array,  r -0,3  X,  testing  wavefront,  M  —  5,  6,  7,  and  8, 


These  results  are  repeated  in  Figure  14  for  the  circular-0  array.  Comparing  Figures  13  and  14  shows 
that  the  circular-0  arrangement  is  generally  inferior  to  the  circular  array  with  the  same  number  of 
elements.  A  heuristic  argument  for  this  result  is  that  the  center  element  measurement  is  the  average  of  all 
the  other  element  measurements  and  does  not  contain  substantially  new  information. 

From  a  collection  of  30  scenarios  that  were  designed  to  represent  realistic  jamming  situations,  three 
cases  were  examined.  The  resulting  eigenvalue  rolloffs  were  compared  to  each  other  and  to  the  results  of 
the  testing  wavefront.  Actual  jammer  locations  and  power  levels  for  these  scenarios  are  given  in 
Appendix  A.  For  the  circular  array  (Figure  15),  the  eigenvalue  rolloffs  for  these  scenarios  are  far  more 
precipitous  than  for  the  testing  wavefront.  They  indicate  that  only  three  or  four  samples  (elements)  are 
useful,  even  though  each  scenario  has  at  least  27  jammers.  Figure  16  shows  the  analogous  results  for  the 
circular-0  array.  The  rolloff  is  slightly  faster  than  for  the  circular  array.  Again,  the  center  element  appears 
to  be  of  dubious  value. 
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344558-14 

Circular  Array,  r  =  0.3 


Figure  14.  Eigenvalue  rolloff for  the  circular-0  array,  r-O.SX,  testing  wavefront,  M  -  5,  6,  7,  and  8. 


344558-15 

r  =  0.3X 


Figure  15.  Eigenvalue  rolloff for  the  circular  array,  r  =  0.3  X,  testing  wavefront,  and  three  jammer  scenarios,  M  —  6. 
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344558-16 

r  =  0.3X 


Figure  16.  Eigenvalue  rollojf  for  the  circular-0  array,  r  =  0.3 X,  testing  wavefront,  and  three  jammer  scenarios, 
M=6. 
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7.  FUTURE  WORK 


When  the  incoming  wavefronts  are  associated  with  significant  multipath,  STAP  processing  is  more 
suitable  and  requires  the  use  of  tapped  delay  lines  at  each  receiver  to  handle  the  correlated  input  [8], 
requiring  replacement  of  the  element  correlation  matrix  with  one  that  includes  the  output  of  the  tapped 
delay  lines.  If  there  are  Af  elements  and  K  taps  per  element,  then  the  correlation  matrix  would  be  an  Mx 
K  hy  M  X  K  matrix.  The  source  or  testing  wavefront  would  need  to  be  modified  to  include  the 
characteristics  of  the  multipath  in  the  wavefronts;  the  eigenvalue  rolloffs  of  this  expanded  matrix  would 
then  be  the  ones  under  investigation. 
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8.  CONCLUSIONS 


For  a  linear  array  of  aperture  Afk,  the  number  of  useful  elements  (samples)  of  practical  utility  is 
limited.  If  the  number  of  elements  is  limited  to  those  spaced  XI2  apart,  they  are  all  useful  and  correspond 
to  2AlX+\  elements.  As  more  and  more  elements  are  added,  the  correlation  matrix  becomes  ill 
conditioned,  implying  small  eigenvalues,  which,  in  turn,  imply  dimensions  in  the  signal  subspace  that 
have  coefficients  too  small  to  be  valuable  as  they  will  be  masked  by  noise.  This  investigation  was 
extended  to  circular  arrays,  with  similar  results.  Clearly,  there  is  a  limit  to  the  number  of  useful  elements 
confined  within  a  fixed  aperture. 
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APPENDIX— EMITTER  AND  JAMMER  DATA  SETS 


%  EMITTER  DATA  SET 

case  Lout 

%  Satellite 

Az  (deg) 

El  (deg) 

Power  (dBW) 

if(cno  =  1) 
SA  =  [  1, 

143.49, 

0.95, 

-161.170; 

2, 

5.83, 

9.17, 

-159.800; 

3, 

106.43, 

64.90, 

-152.330; 

4, 

-15.38, 

67.68, 

-153.830; 

5, 

110.91, 

30.42, 

-158.090; 

6, 

-143.96, 

56.83, 

-154.190; 

7, 

-60.71, 

52.31, 

-157.570]; 

elseif(cno==2) 

%  EMITTER  DATA  SET 

case2.out 

%  Satellite 

Az  (deg) 

El (deg) 

Power  (dBW) 

SA  =  [1, 

61.56, 

15.27, 

-160.920; 

2, 

136.12, 

49.35, 

-159.440; 

3, 

109.83, 

22.89, 

-158.830; 

4, 

-74.70, 

66.43, 

-155.860; 

5, 

-148.32, 

39.24, 

-160.150; 

6, 

0.00, 

12.38, 

-158.670; 

7, 

-13.04, 

34.62, 

-159.880; 

8, 

146.16, 

10.60, 

-159.310]; 

elseif(  cno  = 

3) 

%  EMITTER  DATA  SET 

case3.out 

%  Satellite 

Az  (deg) 

El (deg) 

Power  (dBW) 

SA  =  [1, 

136.42, 

69.88, 

-153.800; 

2, 

-177.46, 

41.47, 

-158.420; 

3, 

^9.71, 

48.92, 

-156.850; 

4, 

-154.15, 

8.14, 

-162.000; 

5, 

13.72, 

8.72, 

-161.380; 

6, 

-15.72, 

0.91, 

-163.230; 

7, 

68.88, 

51.94, 

-156.580; 

8, 

76.82, 

7.24, 

-161.480; 

9, 

-124.97, 

35.40, 

-157.090; 

10, 

-114.35, 

4.76, 

-161.400;]; 

%  11, 

0.0 

0.0 

-130]; 
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%JAMMER  DATA  SET  case  1 

if  (cno  “  1 ) 


%Jamnier 

Az  (deg) 

El (deg) 

Power  (dBW) 

JA-[1, 

56.88, 

54.17, 

-114.154 

2, 

173.15, 

84.09, 

-117.407 

3, 

42.10, 

61.02, 

-136.083 

4, 

53.14, 

55.55, 

-110.831 

5, 

37.69, 

63.86, 

-133.784 

6, 

63.75, 

51.88, 

-124.889 

7, 

49.52, 

57.15, 

-136.222 

8, 

71.28, 

49.96, 

-125.575 

9, 

41.38, 

61.45, 

-135.430 

10, 

53.35, 

55.47, 

-136.617 

11, 

36.85, 

64.41, 

-133.789 

12, 

170.71, 

81.64, 

-117.242 

13, 

174.31, 

85.27, 

-130.050 

14, 

62.51, 

52.20, 

-125.383 

15, 

53.83, 

55.21, 

-115.507 

16, 

99.88, 

48.79, 

-118.332 

17, 

-178.24, 

92.31, 

-122.050 

18, 

177.32, 

88.24, 

-117.136 

19, 

103.32, 

49.15, 

-119.662 

20, 

69.19, 

50.62, 

-113.307 

21, 

71.10, 

49.62, 

-105.939 

22, 

70.14, 

50.22, 

-125.200 

23, 

60.99, 

52.23, 

-106.727 

24, 

172.14, 

82.62, 

-94.359 

25, 

109.48, 

50.29, 

-118.389 

26, 

105.53, 

49.33, 

-119.938 

27, 

104.03, 

49.25, 

-119.507 

28, 

66.36, 

51.12, 

-125.195 

29, 

68.54, 

50.57, 

-125.670 

30, 

56.20, 

54.01, 

-105.094 

31, 

57.60, 

53.80, 

-109.896]; 
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JAMMER  DATA  SET  case2 
elseif(cno  ==  2) 


%Jammer 

Az  (deg) 

El (deg) 

Power  (dBW) 

JA  =  [1, 

-137.93, 

91.43, 

-134.656 

2, 

173.76, 

91.08, 

-122.672 

3, 

-95.60, 

92.69, 

-133.532 

4, 

-119.64, 

92.31, 

-105.596 

5, 

-126.33, 

91.74, 

-134.433 

6, 

175.95, 

90.65, 

-123.497 

7, 

-91.89, 

92.83, 

-137.142 

8, 

-136.62, 

92.44, 

-119.627 

9, 

-130.17, 

92.59, 

-119.741 

10, 

-111.63, 

92.37, 

-106.919 

11, 

-143.27, 

91.07, 

-89.823 

12, 

-131.84, 

92.15, 

-121.165 

13, 

-107.37, 

92.29, 

-136.023 

14, 

-137.33, 

92.04, 

-121.078 

15, 

-101.21, 

92.52, 

-135.178 

16, 

-111.59, 

92.37, 

-132.655 

17, 

-132.37, 

92.72, 

-119.125 

18, 

-170.54, 

90.45, 

-90.893 

19, 

-138.08, 

91.42, 

-134.790 

20, 

171.72, 

90.95, 

-123.440 

21, 

-118.75, 

91.94, 

-135.436 

22, 

-129.76, 

91.51, 

-88.568 

23, 

174.09, 

91.11, 

-122.559 

24, 

122.17, 

91.33, 

-95.775 

25, 

-109.08, 

92.29, 

-134.380 

26, 

172.71, 

91.15, 

-122.376 

27, 

175.66, 

91.10, 

-122.665 

28, 

-116.93, 

92.13, 

-133.236 

29, 

-134.97, 

92.25, 

-120.430 

30, 

-163.70, 

90.52, 

-105.219 

31, 

-121.34, 

91.88, 

-134.551 

32, 

126.37, 

91.48, 

-124.411 

33, 

-102.16, 

92.54, 

-134.209 

34, 

-157.49, 

90.68, 

-104.984 

35, 

-142.02, 

91.55, 

-86.870]; 

35 


JAMMER  DATA  SET  case3 
elseif(cno  “3) 


%Jammer 

Az  (deg) 

El  (deg) 

Power  (dBW) 

JA=[1, 

-76.77, 

92.65, 

-95.774 

2, 

23.26, 

94.46, 

^134.784 

3, 

24.79, 

95.14, 

-119.621 

4, 

29.56, 

94.74, 

-121.159 

5, 

65.81, 

93.48, 

-139.728 

6, 

59.26, 

93.68, 

-134.207 

7, 

44.47, 

94.12, 

-133.237 

8, 

35.04, 

94.26, 

-134.430 

9, 

18.05, 

94.18, 

-89.827 

10, 

-2.45, 

94.10, 

-105.219 

11, 

49.79, 

94.03, 

-106.914 

12, 

40.04, 

94.15, 

-134.553 

13, 

23.41, 

94.46, 

-134.659 

14, 

-72.58, 

93.10, 

-124.408 

15, 

-27.10, 

94.34, 

-123.431 

16, 

42.63, 

94.11, 

-135.441 

17, 

49.82, 

94.04, 

-132.655 

18, 

-26.12, 

94.52, 

-122.364 

19, 

52.32, 

93.87, 

-134.385 

20, 

31.26, 

95.02, 

-119.738 

21, 

-25.06, 

94.48, 

-122.668 

22, 

24.05, 

94.86, 

-121.077 

23 

-22.86, 

94.13, 

-123.494 

24, 

19.34, 

94.45, 

-86.865 

25, 

60.20, 

93.63, 

-135.178 

26, 

41.78, 

94.37, 

-105.597 

27, 

29.06, 

95.21, 

-119.121 

28, 

26.43, 

94.93, 

-120.423 

29, 

-9.30, 

94.04, 

-90.890 

30, 

-24.73, 

94.51, 

-122.550 

31, 

-23.15, 

94.54, 

-122.656 

32, 

31.61, 

94.03, 

-88.574 

33, 

3.78, 

94.17, 

-104.985]; 
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